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ABSTRACT 

We consider models for the generation of the emission line recently discovered in the X-ray afterglow spectrum 
of several bursts, and especially of GRB 99 12 16 observed by Chandra. These observations suggest the presence of 
0.1-1 solar masses of iron in the vicinity of the bursts. We show that there are strong geometrical and kinematical 
constraints on the line emitting material. We discuss several classes of models, favoring one where the line photons 
are produced by reflection off the walls of a wide funnel of semi-opening angle $ « 45° excavated in a young 
(a few months old) supernova remnant made of m 1OM with radius 6 x 10 15 cm, and « 1M Q of iron, providing 
strong support for the SupraNova model. 

Subject headings: gamma rays: bursts — radiation mechanisms: nonthermal — line: formation 



1. INTRODUCTION 

There are now five bursts with evidence for large amounts of 
material around the site of the explosion: four (GRB 970508, 
; Piro et al., 1999; GRB 970828, Yoshida et al., 1999; GRB 
991216, Piro et al., 2000, hereafter P2000; GRB 000214, An- 
tonelli et al., 2000) show an emission feature during the after- 
glow, and one (GRB 990705, Amati et al., 2000) displays an 
edge in absorption during the burst itself. The observations of 
GRB 991216 by Chandra, begun 37 hr after the burst and last- 
ing for about 3 hr, show a 3.49±0.06 keV line significant at 
a ~ 4cr level. If the line is identified with the recombination 
Ka line from H-like iron (6.97 keV), the inferred redshift is 
z = 1.00 ±0.02 consistent with the largest optical absorption 
redshift system (z = 1.02) (Vreeswijk et al., 2000). The line, 
resolved by Chandra, displays a width corresponding to 0.05c 
(P2000). The line and the 2-10 keV continuum fluxes were 
~ 1.6 x 10~ 13 and ~ 2.4 x 10~ 12 erg s" 1 cm" 2 , respectively, im- 
plying (for z = 1.02) an emission rate of iron line photons of 
A/line ~ 4 x 10 52 s" 1 , luminosity Lij ne ~ 4 x 10 44 erg s" 1 and to- 
tal energy £ij ne ~ 3 x 10 49 erg (assuming steady emission for 
40/(1 +z) hr, with H = 75 km s" 1 Mpc" 1 and q = 0.5). If each 
iron atom produces k line photons, the required iron mass is 
M Fe ~ 195 k~ x M Q . Bringing this mass down to ^0.1 M Q 
implies k > 2000 which in turn implies a limit on the recombi- 
nation time 2 . This translates into the condition n e > 10 10 Tj^ 4 
cm" 3 ; the electron temperature T = 10 7 Tj K comes from inter- 
preting the broad bump, marginally detected by Chandra, at 
4.4±0.5 keV as the recombination continuum of H-like iron 
(9.28 keV). Analogously to GRB 950708 (Lazzati et al. 1999), 
the line detection implies the presence of w 0.5M Q of pure iron 
in the vicinity of the burster. 

Observation of the line f Q b s after the burst implies that the line 
emitting material must be within ~ ct b S /(l +z) from the burst 
site: this region must then be compact, contain ~ 0.1M Q only 



in iron, but nevertheless be optically thin to electron scattering, 
such that Comptonization does not broaden the line beyond the 
observed width. We call this the size problem. 

Furthermore, if we interpret the line width as due to the ve- 
locity of a supernova remnant, the limit on the size allows to 
estimate the age of the remnant: for GRB 991216, ~ 15 days. 
At this time (Fig. 1), cobalt nuclei outnumber nickel and iron 
nuclei and the line would be produced mainly by cobalt, not 
iron, at an energy e = 7.5/(1 +z) keV. We call this the kinematic 
problem. 

We discuss in section 2 the constraints implied by the line 
observation; we then present in section 3 three scenarios which 
could produce the observed line. Two of them are found to be 
viable only by adopting ad hoc assumptions or through fine tun- 
ing, while a third scenario, the "wide funnel" scenario, appears 
to be more promising. In section 4 we discuss our results. 

2. PHYSICAL AND GEOMETRICAL CONDITIONS 

The detection of emission features in the afterglow spectra 
of GRBs some hours after the GRB event poses a strong con- 
straint on the location of the line emitting material. If the line 
is detected after f b s from the burst explosion the material must 
be located within a distance R given by 
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where i? is the angle between the line emitting material and the 
line of sight at the GRB site. This limit has important implica- 
tions, because the scattering optical depth is large, 
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(where p = 1, 1.2 and 2 for pure hydrogen, solar composition, 
and no hydrogen, respectively), and, for a remnant radial ve- 

1 Present address: Institute of Astronomy, Madingley Road, Cambridge CB3 0HA, UK; e-mail: lazzati@ast.cam.ac.uk 



2 From the formulae given in Verner & Ferland (1996), 
10 6 < T < 10 8 K. 
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locity of v = 10 9 V9 cm s" 1 the time elapsed from the SN is 
t$N — 12.5(? o i, s /10hr)/[(l +z)(l — cos^Vg] days, implying that 
most 56 Co nuclei (and a fraction of 56 Ni) have not yet decayed 
to 56 Fe (half-lives of 77.3 and 6.08 days, respectively). In Fig. 
1 we show normalized abundances of nickel, cobalt and iron as 
a function of time from the SN explosion. 

To produce the line by frequent recombinations and ioniza- 
tions, a sufficiently high ionizing flux and a large number of H- 
like iron ions, FeXXVI, are required. Efficient use of these ions 
requires that the photoionization optical depth for FeXXVI, 
TFeXXVi IS !• Using photoionization cross section OFexxvi ~ 
1.2 x 10~ 20 cm" 2 , and abundance ratio £ = Mpexxvi/^Fe, we 
have 
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A strict lower limit on the ionizing continuum power comes 
from the line flux itself. For GRB 991216, for which there 
is no rebursting in the X-ray afterglow (contrary to GRB 
970508), assuming an illumination time of 20 hr, we require 
Ziion > 3 x 10 49 erg. Even in the case of an illuminating con- 
tinuum not directly visible by us, we require that the scat- 
tered flux produced by the line emitting material does not af- 
fect the afterglow emission, i.e. it must contribute less than 
10~ 12 erg cm" 2 s" 1 to the observed flux, giving the upper limit 
min(l,T T )£ion < 1.8 x 10 50 erg. 

2.1. Compton broadening and Thomson optical depth 

If the emitting atoms are diluted in an electron cloud, the 
fraction of unscattered flux = [1 — exp(— tt)]/tt, while if the 
scattering cloud surrounds a central source the fraction = 
exp(-Tt). If the scattering electrons have a temperature in the 
~ 10 7 K range (according to the indication that the recombina- 
tion continuum observed in GRB99 1216 may be broad, P2000), 
the average energy shift of line photons scattered only once is 
AE/E « 0.063(/tr/lkeV) 1 / 2 (Pozdnyakov, Sobol & Sunyaev 
1983), already broader than the observed line width. Multiple 
scatterings would further broaden the line, making it increas- 
ingly difficult to detect against the X-ray continuum. If the 
line of GRB 991216 were made by photons escaping unscat- 
tered from a medium with tt > 1 , the problem with the line flux 
would be exacerbated, a factor tj more iron being required. 

An alternative possibility, that the electron temperature <C 
10 7 K, say ~ 10 6 K expected if the ionizing continuum were 
steep and extended to low energies, might conjure with a suit- 
able value of tt to yield, through Compton broadening, the ob- 
served line width. In this case, however, the centroid of the line 
would be redshifted (line photons are mostly backscattered by 
colder electrons) and the recombination continuum narrower. 
Due to poor statistics, the recombination continuum of GRB 
991216 cannot place a firm constraint on the electron tempera- 
ture, but the models in Section 3 are largely independent of the 
presence and width of the recombination emission continuum. 

3. THE MODELS 

We can roughly divide the models into transmission and re- 
flection models. If we see the line in reflection, line photons 
come from the layer with rp e xxvi = several 3 . For an iron abun- 
dance greater than the solar value, in this layer tj < 1 . In gen- 
eral, transmission models require tj < 1 and Tp e xxvi ~ 1, while 

3 What we call FeXXVI might well be CoXXVII or NiXXVIII, see section 2. 



reflection models require TFexxvi > tt. Each iron atom must 
produce ~ 2000 iron photons, and this requires an electron den- 
sity larger than 10 10 cm" 3 . 

The large equivalent widths of the X-ray line in GRB 991216 
(EW = 0.5 keV) and, especially, GRB 970508 (EW ~ 1 keV), 
GRB 970828 (EW ~ 3 keV) and GRB 000214 (EW ~ 2 keV) 
favor models in which the line originates in reflection, not trans- 
mission. In reflection, in fact, the ionizing flux is always effi- 
ciently reprocessed in line photons in the rp e ~ 1 layer, while 
in transmission this happens only for a particular tuning of the 
FeXXVI and the free electron densities. 

Three alternative geometries (Fig. 2) producing lines in re- 
flection are described below. 

In reflection models, we can derive the photon line lumi- 
nosity by estimating the volume V em effectively contributing 
to the line emission, and assuming a given iron mass. If the 
layer contributing to the emission has tj ~ 1 (in order to avoid 
excessive Compton broadening of the line), and in this layer 
rpexxvi few (to efficiently absorb the continuum), we have 
Vem = S/(<7Tn e ), where S is the emitting surface. The line emis- 
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where we assumed that the total volume is V = SAR (slab or 
shell geometry). Eq. Q shows that the total iron mass must be a 
sizable fraction of a solar mass to obtain the observed line flux 
(4 x 10 52 s" 1 ), but less than the largest observed values, 0.9M Q 
for SN 1998bw (Type Ic, Sollerman et al., 2000) and SN1991T, 
(Type la, Filippenko et al., 1992). Eq. ||all by itself establishes 
that the line emitting material must be a SNR: no other known 
astrophysical object contains so much iron, not even the largest 
known stars (though these fail by a factor of a few only). 

Eq. [|has another important implication: it fixes M/R, and 
since R is fixed (for a given geometry, see below) by the time- 
delay, the cloud density is also fixed, to within an order of mag- 
nitude. Solutions with densities much exceeding (or much be- 
low) « 10 10 cm" 3 are excluded. 



3.1. The wide funnel 

If the GRB explodes within the young remnant of a super- 
nova, we must consider two possibilities: a plerionic and a shell 
remnant. In the first case, we consider a wide funnel excavated 
in a young plerionic remnant. This can solve the size prob- 
lem, since it extends to large radii but, at the same time (Eq. [I]), 
can maintain the time-delay contained because it is naturally 
built close to the polar axis (assumed close to the line of sight). 
The geometry is sketched in Fig. 2a. Fixing the line photons 
rate (Eq. yields R = 6x 10 15 cm, and thus an opening angle 
•d = 48°, to fit the time-delay. For any reasonable SN composi- 
tion, these parameters imply « e > 10 10 cm" 3 . The funnel walls 
are probably not straight (like in a cone geometry), but curved 
instead, like the surface of the coffee in a cup when it is spun 
up, and can be illuminated by the ionizing flux from a central 
source. Assuming a cone geometry for simplicity, we rewrite 
Equation Ej as 
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where in this case R (the radius of the filled remnant) has been 
used instead of AR and tan$ takes into account the geometry 
of the funnel. This is a lower limit, since a parabolic geometry 
has a larger funnel surface and we neglected the (likely) density 
stratification inside the remnant. 

Since the funnel's normal is not parallel to the incident pho- 
ton's momentum, radiation pressure exerts a force parallel to 
the surface on the layer tj < 1 (which reflects isotropically pho- 
tons arriving from a specific direction). Calling ip the angle 
between the funnel's normal and the incident photon's arrival 
direction the absorbed fluence E ion accelerates the funnel layer 
to a speed Vf = smip(2Ei on /Mi aye[ ) 1 / 2 ~ 10 4 £^ 2 50 sin^ km s -1 

when R = 6 x 10 15 cm. Thus, ablation by radiation pressure 
propels the reflecting layer to velocities like those seen in 
GRB991216byP2000. 

The funnel model can thus explain the observed line broad- 
ening even in a relatively old SNR. For a remnant with v e j = 
3000 km s" 1 , a typical value for type II SNe, we find a remnant 
age fsNR ~ 230 days, enough to turn most 56 Ni and 56 Co into 
56 Fe (Fig. 1). The emitting iron is then boosted to higher veloc- 
ity by the ionizing flux, circumventing the kinematical problem. 

3.2. Back of the remnant 

We consider now a shell remnant, characterized by a wide 
opening, so that a substantial fraction of its inner face is visi- 
ble (Fig. 2b), and illuminated by burst and afterglow and, more 
likely, by the shock due to the fireball impact (Bottcher 2000). 
Line photons are produced in the layer Tf e xxvi ~ few. As be- 
fore, the remnant may have been expanding slowly, until accel- 
erated by burst photons. Since we see the remnant's internal 
part, the size limit discussed in Section 2 is tighter, and we de- 
tect mainly redshifted photons. For a remnant distance R = 10 15 
cm and velocity 5,000 km s" 1 , the age is ~ 20 days, close to the 
56 Co peak (Fig. 1). However, the 7.5 keV 56 Co line would be 
observed with a substantial redshift (and therefore mimic a 56 Fe 
6.97 keV line), if the burst photons increased the expanding ve- 
locity to ~ 20, 000 km s" 1 (implying that the burst emitted a few 
times 10 52 erg sideways, for a 10 M Q remnant). In this model 
the observed line width of 15,000 km s" 1 requires that only a 
limited range of projections of the radial expansion velocity 
vectors are observed. This solution is not convincing since it 
involves fine tuning between the energy produced sideways and 
the mass of the remnant, in order to have the correct expansion 
velocity, even though it explains naturally the different centroid 
energy of the line observed in GRB 970508: since the line is 
redshifted by shell expansion, we expect a distribution of ob- 
served energies as a consequence of different velocities of the 
remnant. 

3.3. Back-illuminated equatorial material 

We now explore the possibility of a simultaneous GRB- 
"supernova" explosion. Consider a scenario in which the GRB 
ejects and accelerates a small amount of matter in a collimated 
cone, while a large amount of matter is instead ejected, at sub- 
relativistic speeds, along the progenitor's equator. For illus- 
tration, assume equal amounts of energy, 10 52 erg, in both di- 
rections, and 1 Mq expelled along the equator. Then v eq /c = 

Q.\A%e\!^{M/Mq)~ 1 I 2 . Massive star progenitors are always 
surrounded by dense material from strong winds with mass 



loss rates m w = 10 m w ._5 and velocity v w = 10 7 v w .7. This wind 
scatters back a fraction of the photons produced by the bursts 
and its afterglow (Thompson & Madau 2000), ,nd, in return it 
is also accelerated until it reaches relativistic velocities, when 
scattering efficiency decreases. Thus assume that each elec- 
tron (and the associated proton) scatters photons only until it 
reaches Y = 2, i.e. until it has scattered a total energy of m p c 2 
(i.e. ^2,000 photons of 0.5 MeV each). In this case the scat- 
tered luminosity L scatt is constant, since there are an equal num- 
ber of electrons in a shell of constant width AR (for a density 
profile oc R~ 2 ). This luminosity is of the order: 

t 2 m w i o 1 rv45 m w.-5 -1 ,,, 

iscatt ~ m p c z = 1.8 x 10* 3 ergs . (6) 

TOpV w /C V W J 

Scattered photons illuminate the expanding equatorial matter 
after a time 2R / c, producing X-ray line emission mostly due to 
Ka emission from H-like 56 Ni at 8.1 keV. Even if the large ex- 
panding velocity makes the transverse Doppler shift important, 
this requires a fine tuned velocity to redshift the 8.1 keV 56 Ni 
line into 6.97 keV, for which reason we consider this scenario 
unlikely. 

4. CONCLUSIONS 

X-ray line emission features from GRBs' afterglows give us 
information on the nature of their progenitors by imposing con- 
straints on models, the most severe being how to arrange a large 
amount of line emitting material around the GRB site, but with- 
out large Thomson scattering opacities. A further limit comes 
from the line width of Chandra's observation in GRB 991216. 
These observations require a large iron mass (Eq. Q), such as is 
seen only in SNe. 

We have presented three distinct models: two (the "back of 
the remnant" and "back illuminated equatorial material" mod- 
els) require that the observed emission line is produced by red- 
shifted cobalt or nickel (instead of iron), and need fine tuning 
and/or ad hoc assumptions. The "wide funnel" model, instead, 
solves the size problem, and the acceleration of the line emit- 
ting material by grazing incident photons solves the kinematic 
problem, allowing time for cobalt to decay to iron. It implies 
that the GRB progenitors are massive stars exploded as super- 
novae months before the burst, contaminating the circum-burst 
environment with iron rich material. This two-step process, 
and the time-delay between the two steps, are predicted by the 
SupraNova scenario (Vietri & Stella, 1998). Alternative expla- 
nations invoke extremely iron enriched massive winds (Weth et 
al. 2000). 

Future observations may address issues such as the presence 
of cobalt and nickel lines in the spectrum, the age of the rem- 
nant, the line profile, the geometry and kinematics of the emit- 
ting region, the presence of more than one line, the ionization 
state, the time evolution of the line and edges, all indicative of 
the characteristics of the illuminator. Together, this information 
may determine whether the line emitting material originates in 
a normal, albeit rare, supernova remnant or if a more exotic 
explosive phenomenon is required. 

We thank Frits Paerels and Martin J. Rees for useful discus- 
sions. 
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Fig. 1 . — The normalized abundances of 56 Ni, 56 Co and 56 Fe as a function of time from the creation of 56 Ni. The neutral, He-like, 
H-like and recombination edge energies are (in keV): 56 Ni: 7.478, 7.806, 8.102, 10.775; 56 Co: 6.930, 7.242, 7.526, 10.012; 56 Fe: 
6.404,6.701,6.973,9.278. 




Fig. 2. — Cartoon of the three models discussed in this paper. In panel a) a funnel excavated in a supernova remnant produces the 
iron line in reflection, by material accelerated to ^10,000 km s" 1 by the strong radiation pressure. In panel b) we see the internal 
parts of a young and receding supernova remnant, in panel c) some equatorial material exploding simultaneously with the burst is 
illuminated by burst and afterglow photons scattered by the pre-burst stellar wind. 



